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Abstract 

A one parameter, model confined-gluon propagator is employed in a phe- 
nomenological application of the Dyson-Schwinger and Bethe-Salpeter equa- 
tions to the calculation of a range of n- and p-meson observables. Good 
agreement is obtained with the data. The calculated quark propagator does 
not have a singularity on the real-p 2 axis. A mass formula for the pion, in- 
volving only the vacuum, dressed quark propagator, is presented and shown 
to provide an accurate estimate of the mass obtained via a direct solution of 
the Bethe-Salpeter equation. 
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I. INTRODUCTION 



The Dyson-Schwinger equations (DSEs) provide a useful, semi-phenomenological tool for 
the study of QCD. These coupled integral equations relate the n-point (Schwinger) functions 
of QCD to each other. They provide a nonperturbative, Poincare invariant framework 
that enables one to correlate hadronic observables through the properties of the Schwinger 
functions of the elementary excitations in QCD; i.e., the Schwinger functions of quarks and 
gluons. (Quark and gluon propagators (2-point functions) are examples of such Schwinger 
functions.) This makes it particularly suitable for addressing questions such as confinement 
and dynamical chiral symmetry breaking and also hadronic spectroscopy and interactions. 
This approach is reviewed in Ref. |] and has recently been applied to the study of ir-n 
scattering 0, the electromagnetic pion form factor M, p-u mixing M and the anomalous 
7*7T° — ► 7-transition form factor. || 

It is possible to obtain information about such Schwinger functions via a numerical sim- 
ulation of a lattice-QCD action. |6|-Q However, in addition to the usual problems associated 
with identifying and establishing the existence of the continuum limit, and recovering the 
global symmetries of QCD, this also requires gauge fixing on the spacetime lattice. Gauge 
fixing eliminates a number of gauge-equivalent gauge-field configurations, thereby leading 
to poorer statistics. It does not eliminate all such configurations, however. One is left with 
Gribov copies; i.e., gauge configurations in the gauge- fixed ensemble that are not distinct 
but are related by topologically nontrivial gauge transformations. [0HTT| This entails an over- 
counting problem in the evaluation of gauge-fixed correlation functions. Present studies are 
encouraging, having established that this approach to the calculation of gauge-fixed QCD 
Schwinger functions is feasible. |7| However, the problems identified above entail that they 
are currently qualitatively and quantitatively unreliable. 

Presently the most reliable estimates of the behaviour of quark and gluon Schwinger 
functions are obtained in DSE studies. The DSEs are a tower of coupled equations and a 
solution is only tractable if this tower is truncated. Truncation procedures that preserve 
the global symmetries of QCD are easy to construct and implement. This has not yet been 
accomplished for the local symmetry in QCD, however, progress is being made following 
the realisation of the importance of the nonperturbative structure of the fermion-gauge- 
boson vertex. [^|T2|-|iT)H This introduces an uncertainty in the infrared; i.e., for k 2 < 1 — 
2 GeV 2 . However, this uncertainty is merely quantitative. There is general agreement on 
the qualitative features of the quark and gluon 2-point Schwinger functions; i.e., 1) that the 
gluon 2-point function is significantly enhanced at small spacelike-fc 2 [ji],pl7)-|T9|j and that this 
entails an enhancement of the momentum-dependent quark mass-function [|I|,f2l^-|24| ; and 2) 
that for k 2 > 1 — 2 GeV 2 the two-loop, renormalisation group improved, perturbative results 
are quantitatively reliable. 

Some phenomenological DSE studies have employed a parametrisation of the 2-point 
quark Schwinger function based on these results; for example, Refs. [0-0. Such studies are 
phenomenologically efficacious. However, they involve the addition of new parameters when 
applied to systems involving other than u and d quarks. 

The introduction of new parameters is unnecessary when the propagator for a quark of 
a given flavour is obtained directly from a quark DSE whose kernel is determined by the 
2-point gluon Schwinger function and the quark-gluon vertex. This procedure correlates the 
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propagators for quarks of different flavours via the parameters in the gluon 2-point function. 



There have been studies that employ this approach; for example Refs. pT|j25|j26[| . However, 
it is computationally more intensive and the studies therefore addressed the calculation of 
a smaller class of observables. The present study is a first step in extending this latter 
approach. 

Herein we employ a one parameter model gluon propagator (gluon 2-point Schwinger 
function), motivated by the results of Refs. O-O], in a calculation of a range of tt- and 
p-meson observables. The one parameter is a mass scale that can be interpreted as marking 
the transition between the perturbative and nonperturbative domains. This model gluon 
propagator provides the kernel for a quark DSE, which is solved to obtain the quark prop- 
agator (quark 2-point Schwinger function) for real-p 2 G (—00,00). These two Schwinger 
functions provide the kernel of the meson Bethe-Salpeter equation (BSE), whose solution 
yields the meson mass and Bethe-Salpeter amplitude, which is a necessary element in the 
calculation of decay constants and scattering lengths, for example. The single mass param- 
eter determines all of these Schwinger functions and is varied to obtain a good fit to a range 
of calculated n observables. This illustrates the utility and economy of the approach. 

In studying the pion BSE we derive a mass formula for the pion, which involves only 
the vacuum, dressed quark propagator, valid to all orders in rriR, the renormalised current 
quark mass. Our numerical studies show that this formula provides an excellent estimate of 
the mass that is obtained by actually solving the BSE. 

The model gluon propagator is discussed in Sec. and the quark DSE in Sec. |T|. The 
pion mass formula is presented in Sec. [TV|- Our numerical results are discussed in Sec. [V] 
and we summarise and conclude in Sec. |VT[ 



II. MODEL GLUON PROPAGATOR 



In Euclidean metric [27] the Landau gauge gluon propagator is 



g2D ^ k) = - ^ P[1 + H(P)] (1) 

where U(k 2 ) is the gluon vacuum polarisation. Setting Zf h = Z§ h , where Zf h is the renor- 
malisation constant for the ghost-gluon vertex and Z§ that for the ghost wave function, 
then 

satisfies the same renormalisation group equation as the QCD running coupling constant, 
a(k 2 ), |28| and hence 

(^ W ) H =(V-^)^. (3) 

This is sometimes described as the "Abelian approximation" because it entails the QED-like 
Ward identity Z\ = Z2, where Z\ is the quark-gluon vertex renormalisation constant and 
Z2 is the quark wavefunction renormalisation constant. 
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The two-loop renormalisation group expression for the running coupling constant only 
receives small (~ 10%) corrections from higher orders for spacelike-A; 2 > 1 GeV 2 and hence 
can be said to provide an accurate representation on this domain. For k 2 < 1 GeV 2 , however, 
a(k 2 ) is not known and can only be calculated nonperturbatively. The current status of such 
studies is summarised in Ref. |[] and, as remarked in Sec. [I], gluon-DSE studies agree on the 
qualitative behaviour of a(k 2 ) at small-A; 2 . Present phenomenological quark-DSE studies 
rely on an Ansatz for a(k 2 < 1 GeV 2 ) motivated by these gluon-DSE studies. 

Herein we consider a parametrisation suggested by the Landau gauge studies of Ref. pjl . 



which revealed a strong enhancement in the gluon propagator at small spacelike- k (< 
1 GeV 2 ) that could be described by an integrable singularity. We employ the one parameter 
form: 



A(k 2 ) =An 2 d 



A7r 2 mi5\k) + 



1 _ e (-fc 2 /[4m?D 
k 2 



(4) 



where d = 12/(33 — 2Nf), with Nf = 3 the number of light flavours. The first term in 
Eq. (^) provides an integrable, infrared singularity p0| , which generates long-range effects 
associated with confinement, and the second ensures that the propagator has the correct 
large spacelike-/c 2 behaviour, up to ln[fc 2 ]-corrections. A form similar to this has been used 



by other authors [pl|-p4]| with 1-loop logarithmic corrections included in the second term. 
We neglect these terms as a simple expedient to ensure that our gluon propagator does not 
have a Lehmann representation and may therefore be interpreted as describing a confined 
particle; i.e., an elementary field with which there is no associated asymptotic state. |l|j29|| 
Since ours is a model gluon propagator there is no reason why the coefficients of the 
two terms in Eq. @ should be related in the particular fashion we have chosen. However, 
consider 



A fx 2 



d A k 



e ik - x A(k 2 ) = d 



m 2 + ±-e- x2m < 
x 2 



(5) 



It is clear from this that with our choice of the ratio of these coefficients the effects of <5 4 (/c) 
in Eq. (0) are completely cancelled at small x 2 ; i.e., 



A(x* 



i 2 x 2 <i d 



3 + 



(6) 



which is the form expected from QCD (again neglecting logarithmic-corrections). 

One can therefore interpret m t as the mass scale in our model that marks the transition 
from the perturbative to the nonperturbative regime. Herein m t is varied to provide a best 
fit to a range of calculated pion observables. [See Eq. (|52|) and the associated discussion.] 



In Euclidean metric 



III. QUARK SELF ENERGY 

the DSE for the quark propagator is 
S~ l {j)) = Z 2 (vy -p + mo) + E'(p), 



(7) 



where 
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V{p) = z if-0jil9 2 d, u (p - k)^S(k)r„(p, k) , (8) 

with T^(p, k) the quark-gluon vertex, is the regularised self energy, which can be decomposed 

as 

E'(p)=H-p(^V)-l)+B'b 2 ) • (9) 

The inverse of the renormalised quark propagator is 

S- 1 (p)=i 1 -p + i:(p)=i 1 -pA(p 2 ) + B(p 2 ) . (10) 

Herein the prime denotes regularised quantities and unprimed quantities are fully renor- 
malised. 

We employ a subtractive renormalisation scheme, requiring that, at a spacelike renor- 
malisation point, fj 2 , 

5 _1 (p)| p 2 =M 2 =inf-p + m R , (11) 

with m R the renormalised current quark mass. 

In this scheme, the wavefunction and mass renormalisation constants are given by 

Z 2 = 2- A\fi 2 ,A 2 ) and m R = Z 2 m (A 2 ) + B'(fi 2 , A 2 ) , (12) 

respectively, and the renormalised self energies are therefore obtained from 

A(p 2 , v 2 ) = 1 + A'(p 2 , A 2 ) - A'Qm 2 , A 2 ) , (13) 
B(p 2 , /i 2 ) = m R ^ 2 ) + B'(p 2 , A 2 ) - 5 V, A 2 ) . (14) 

In this scheme, A(fx 2 ) = 1 and -B(/i 2 ) = m R (fi 2 ). (In the following we often write m R (n 2 ) as 
simply m R , in which case the \j? dependence is implicit.) 
The renormalised axial-vector Ward identity is 

(P - QU r>, q) = S-\p) 75 + 75 S'\q) - 2m R T 5 (p, q) . (15) 

The composite operators and T 5 are renormalised such that, at p 2 = fi 2 = q 2 , 
r 5(P> 0) = 7^75 and Y 5 (p, q) = 75. 

The chiral limit is identified as the limit in which the renormalised axial-vector current 
is conserved; i.e, with the limit m R (fi 2 ) — > 0. 

A. Analysis of the large-p 2 behaviour of the quark propagator 

At large spacelike-A; 2 and p 2 one may replace the gluon propagator and the quark-gluon 
vertex by their asymptotic forms: 

A(A; 2 )^-i and r>, k) -> 7 „ . (16) 

In this limit A(p 2 ) = 1 and B(p 2 ) is the solution of 
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B{x) = Z 2 m + - f dyy (- 9{x - y) + -9{y - x)) _f^. (17) 
4 Jo \x y J y + B\y) 

where x = p 2 , y = k 2 and A = 4 Z\ d. 

For x such that B(x) 2 <C x; i.e., for x > /i 2 , this integral equation is equivalent to the 
differential equation 

^ 2 ^ b(i) K b(iHo ' (i8) 

subject to the boundary condition 

B(fi 2 ) = m R (19) 

or 

4-[xB(x)U =Z 2 m . (20) 




\dx 

Under the change of variables x = fi 2 exp(2z), Eq. ( pp becomes 

B(z) + 2S(«) + XB(z) = , (21) 

which is the equation of motion for a damped harmonic oscillator. One has critical damping 
for A = Ac = 1 and this yields the critical coupling for dynamical chiral symmetry breaking; 
i.e., in the absence of the first term in Eq. @, the model would still exhibit dynamical 
chiral symmetry breaking for A > 1. This behaviour has been observed in QED JHJ and 
phenomenological models of QCD without an infrared- singular model gluon propagator. 



The solution of Eq. (|18D consistent with Eq. (|T9|) is 

B(z) = k e~ z cos (z V\ - 1 + 0) , (22) 

with 

K cos 4> = tur . (23) 

In the chiral limit m R = and hence <p = ir/2. In general k is only determined in a complete 
solution of the integral equation. 

The boundary conditions in Eqs. ([19]) and (^TJ) are equivalent: a given value of m# 
entails a given value of Z 2 mo and vice-versa. In fact, for finite A, = generally entails 
■^2^0 7^ 0. It follows from Eq. ([22]) , however, that for any finite value of m# 



lim Z 2 (/i 2 ,A 2 )m (A 2 ) = . (24) 

A 2 ^oo 

Equation (p2|) indicates that the renormalised mass function will exhibit damped oscil- 
lations about zero for p 2 > /1 2 , a feature we observed in our numerical solutions, which were 
well described by Eq. (|22"|) on p 2 e [/i 2 ,A 2 ]. With the exception of Ref. [34j], other DSE 
studies implicitly use /i = A and hence the oscillations are not observed. The oscillations 
were observed in Ref. |[34j| , which addresses in detail the nonperturbative renormalisation of 
the fermion DSE in QED. 
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B. Additional remarks 



The "Abelian approximation" entails that Z\ = Z 2 in Eqs. (|7|) and (H). We make this 
identification hereafter. 

In the numerical studies described below we employed the rainbow approximation: 

r M (p,*) = 7M- (25) 

This is a quantitatively reliable approximation in Landau gauge. (This is not the case in 
other gauges). For example, in studies of the critical coupling for dynamical chiral symmetry 
breaking, a comparison of the results obtained using this approximation |30| with those 



obtained using more realistic vertex Ansatze [|T|, p~2|JT~5| , [T6|j35H shows this approximation to 
be accurate to 5%. The improvements to this approximation are qualitatively important 
Jl, 12, 15, 16, 35], being crucial to the restoration of multiplicative renormalisability and gauge 



covariance. However, herein a quantitatively reliable calculation scheme is sufficient and this 
is provided by Eq. (j25"D in Landau gauge. 

IV. A PION MASS FORMULA 

The unrenormalised BSE for the pion in generalised-ladder approximation is, with un- 
renormalised n-point functions beyond denoted by ~, 

I\r(p; P)+J 44l I 9 2 D^( P - q) 7m S(q + \P) f P) S(q - \P) 7 „ = , (26) 

where P = p\ + p 2 is the total momentum and p = (pi — P2V2 the relative momentum of 
the q-q pair. 

For the pion it is a good approximation [26,[36| to write 



T n (p;P) = l5 F(p 2 ,P 2 ) , (27) 

in the sense that T n (p; P) is a general pseudoscalar 4x4 matrix and the right-hand-side is, 
pointwise, a good approximation to it and the inclusion of the other allowed Dirac amplitudes 
alters the mass eigenvalue by < 1 %. With this approximation Eq. ( |2"ED becomes [C 2 (i?) = 
{N 2 C - l)/(2N c ) = 4/3 for N c = 3] 



8N c F(p 2 ,P 2 ) = 3C 2 (R) J «A(p-q)H(q;P) (28) 
with 

H(p; P) = 8N C (p + ■ p_a+dy + &p s ) F(p 2 , P 2 ) , (29) 
where we have defined p± — p ± P/2, 

50 = - Ml and ^ = . Ml . (30) 

v p 2 ± A{pif + B{p 2 ± y 3 piA{piY + B{piY y 1 
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Equation (p8| ) is a convolution in four dimensions and can be rewritten as 

= 8N c F P (x) -C 2 (R)2,k(x)H P (x) (31) 

with Hp(x) the Fourier transform, with respect to p, of H(p; P). 

Multiplying the right-hand-side of Eq. ([31] ) by (Fp(— x)/[3C 2 (R)A(x)]) one can construct 

^-/-^(s^^^-^H- (32) 

In the auxiliary-field bosonisation of the Global Colour- symmetry Model JI],|T7|] the effective 
action contains the term 

| d 4 x A ^(x) IL^x - y) n%y) , (33) 

with it 1 (x) a local field variable identified with the pion field. One sees from this that U n (P) 
plays the role of the inverse propagator for the composite pion field. Further, at the solution 
of the BSE, P 2 = — m 2 , Eq. (|3T|) is satisfied and hence 

IT^P 2 = -ml) = . (34) 

It has been shown [3S| that for mo = the unrenormalised BSE has a massless, P 2 = 0, 
solution with 

F P (x) = F P =o(x) = B mo=0 (x) , (35) 

which is the manifestation of Goldstone's theorem in the DSE approach. Using this as an 
approximation for P 2 = —m 2 ^ 0, via the unrenormalised DSE: 

B mo=0 (x) = 3C 2 (R) A(x) a^=\x) , (36) 

one obtains 

U n (P) w J d'xB^oix) (8N c a£ 0=0 (x) - H P (x)) = U n (P) . (37) 

This is manifestly invariant under renormalisation and hence one may write 

n*(P) = J d 4 xB mR=0 (x){8N c af R=0 (x)-Hp(x;m R )) , (38) 

with every quantity on the right-hand-side renormalised (as and H have the same form but 
with unrenormalised quantities replaced by renormalised ones) and evaluated with mp ^ 
unless otherwise specified. 

As remarked above, U(P 2 = —m 2 ) = at the solution of the BSE. Equation (5S) 



therefore allows one to obtain a simple pion mass formula derived from the generalised- 
ladder approximation to the BSE and expressed solely in terms of the massless and massive 
renormalised, vacuum, dressed quark propagators. 

For the pion (because m 2 ~ 0) it is a good approximation to write 
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fl n (P) ^ 11,(0) + P 2 iV 2 (39) 

where 
N c r« 



8tt 



/■A 2 

ds s B mR=0 (s) 2 (a v - 2 [g s g' s + soyoy] - s 



II ( I \2 



s 2 



ii i i \2 

ff^y - (cr v ) 



with the primes denoting differentiation with respect to s = p 2 and oy and as evaluated at 
m^. This is just the conventional, generalised-ladder approximation Bethe-Salpeter ampli- 
tude normalisation constant, calculated neglecting small (~ 2%) 0(m 2 ) corrections. 

We note that if A(p 2 ) = 1, N n = f n . In general, the approximation N w f n is accurate 
to within 10% and the difference is a measure of the error introduced by the approximation 
of Eq. (£7]). § (Also see Table. § 



Equation (|3^) yields the explicit pion mass formula |39] 

ml N 2 = A f A2 ds s ( BmR Ms) *r = \s) - B mR=0 (s) . (41) 

One notes immediately that, for a given value of m^, m 2 — > constant < oo as N c — > oo 
and that, for arbitrary Af c , m 2 — > as — > 0. Further, if the DSE is solved with a quark- 
gluon vertex that ensures multiplicative renormalisability then m 2 is a renormalisation point 
invariant and the result is independent of the cutoff A 2 . The integral on the right-hand-side 
of Eq.fl4"T|) is convergent in the limit A 2 — ► oo. 

From Eq. ([0]) one can recover what is sometimes called the Gell-Mann-Oakes-Renner 
relation in the form: 



rnlti = -m£(qq)C, (42) 



where 



N M 2 

-(^)vac = ^ I dssa™«=\s), (43) 

which is the customary definition of the vacuum condensate. However, in terms of the 
nonperturbatively dressed quark propagator, equality between the integrands requires the 
following ad hoc and mutually incompatible "approximations": Vs, 

B mR=0 (s) « B mR ^ (s) ; (44a) 
°s R= \s) » a^°(s) ; (44b) 

B mR ^ (s) ^m R + B mR=0 (s) , (44c) 

which yields Eq. (|^) when one makes the additional approximation N n « /„., discussed 
above. That these are bad "approximations" is clear; for example, Eq. ( |44a| ) has the effect 
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of replacing a convergence factor in the integrand by unity and it is incompatible with 
Eq. ([He]). As elucidated in Ref. [fi0|1 , Eq. (f42D can only be obtained if the (renormalised) 
current quark mass is treated strictly as a perturbation. The inadequacy of Eqs. ( f4"2| ) and 
((43|) is only exposed by a careful treatment of the Dyson-Schwinger and Bethe-Salpeter 
equations. 

We emphasise that Eq. (f4lD is completely consistent with the general arguments of 
Ref. | j41{| . It is derived from the generalised ladder BSE and measures the expectation value 
of the explicit chiral symmetry breaking term in the pion state under the approximation 
that Eq. ( |3"5] ) is valid for P 2 ^ 0, which is why the right- hand- side involves only vacuum 
quantities: massless and massive, renormalised, vacuum, dressed quark propagators. 

We demonstrate below that Eq. ( f4~T|) provides an extremely accurate estimate of the pion 
mass obtained by solving the pion BSE in generalised-ladder approximation. (See Eq. ( |53|) 
and Table |.) 

A. Solving the pion Bethe-Salpeter equation. 

In our numerical studies we are interested in the subtractively renormalised Bethe- 
Salpeter amplitude, F(p; P). This is defined in terms of the regularised amplitude F'(p; P) 
via 

F{ P] P) = F\p-P)-F'{^P) , (45) 
which, in generalised ladder approximation, is obtained as the solution of 

F'(p; P) = Z 2 3 C 2 (R) J A d 4 qA(p - q) (q + ■ q_ a v a v + o%Og ) F(q; P) . (46) 

It is clear that all corrections to free-field behaviour vanish at the renormalisation point; 
i.e., F(p;P)\ p2=ii2 = 0. 

Upon comparison with the DSE for B(p 2 ) in Sec. |Tj, ^ i s clear that in the chiral limit 
(rriR = 0) one has 

F(p;P)=B rnR=0 (p) ; (47) 
i.e., that Goldstone's theorem is manifest. [[38 



One may solve Eq. (^) numerically by introducing an eigenvalue, A(P 2 ), on the right- 
hand-side. This yields an equation that has a solution at every value of P 2 . The equation 
can then be solved repeatedly until that P 2 is found for which A(P 2 ) = 1. 

The eigenvalue and eigenvector are determined by employing the Tschebyshev decom- 
position 



F(p; P) = J2 F *(P\ P 2 ) ^(cos/3) (48) 
i=i 

and solving for the Tschebyshev moments of F(p; P), which are obtained via 

F^p 2 , P 2 ) = 1 [* dp sin 2 p Ui(cosP) F(p, P) . (49) 

J 
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In practice we only keep the lowest moment Fo(p 2 , P 2 ); neglecting the coupling to the higher 
moments. This is a very good approximation for the pion. [P6|| 

For an on-shell pion P 2 = —m 2 and hence the right-hand-side of Eq. (|46| ) samples the 
quark propagator at complex values of its argument. To avoid solving the quark DSE off the 
real-p 2 axis we expanded (q + ■ g_ <JyOy + to 0(P 2 ) and solved the resulting equation, 

which involves derivatives of the propagator at real-p 2 > 0. 



V. NUMERICAL RESULTS AND PHENOMENOLOGY 

We have two parameters: the mass scale m t in the gluon propagator, which marks the 
transition point between the perturbative and nonperturbative domains, Eq. (§); and m R , 
the renormalised current quark mass. We varied these parameters in order to obtain the 
best x 2- fit to the pion observables: [calculated using Eq. (|41|)1, the weak pion decay 
constant M 

^ = An 2 Jo dss ^ F °( s ' P ^ [ av(Js + ^ s ( a 'v a s - vvv's)\ , (50) 

r n and the n-n scattering lengths: a°, a 2 ,, a\, a®, expressions for which are given in Ref. 0. 

At each pair of parameter values the quark DSE was solved numerically with /i = 
48 fm _1 = 9.47 GeV, which is large enough to be in the purely perturbative domain, and 
A = 2 18 fm _1 ~ 5461/i. The results were almost independent of the cutoff; doubling it 
leading only to a 3% change in f n , for example. Our results would have been completely 
independent of A if we had employed a vertex that preserves multiplicative renormalisability. 
This observation provides a quantitative measure of the violation of multiplicative renor- 
malisability when the rainbow approximation is used in Landau gauge. It is significantly 
worse in other gauges. As remarked above, rainbow approximation entails a loss of gauge 
covariance. Our experience suggests that our results would change by no more than 10% 
if we had used a dressed fermion-gauge-boson vertex that ensured gauge covariance of the 
fermion propagator. flI|, |15|L6 ,35 



The formulae for the observables were then evaluated using the solution obtained and 
the approximation that Eq. (fl7| ) is valid for uir ^ 0. After obtaining the optimal values 
of the parameters we recalculated the observables using the pion Bethe-Salpeter amplitude 
calculated as described in Sec. [IV A| . We found numerically that 

F^-^P^^B^p 2 ) . (51) 

The best x 2- fit was obtained with 

m t = 0.69 GeV and m R = 1.1 MeV . (52) 

We also carried out an extended x 2_ fit where the ratio of the coefficients of the two terms 
in Eq. (^) was allowed to vary. In this case the best x 2 was obtained with the value of m t 
in Eq. (|52"D and a ratio that agreed with that in Eq. (J4|) to within 2%. The data therefore 
requires both terms in the propagator and the cancellation of long-range effects described 
in Eq. (H). 
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The observables calculated with these parameter values are presented in Table |. One 
observes immediately that our one parameter model for the gluon propagator provides a good 
description of low energy pion observables. This improves upon the results of Refs. 0-0 , 
in which the quark propagator was parametrised and illustrates the connection, suggested 
in these articles, that may be made between hadronic observables and the quark-quark 
interaction. 

We have made a direct comparison on the spacelike-p 2 axis of the numerical solutions 
for ay and a$ obtained herein with the parametrised forms used in Ref. |J. The agreement 
in form and magnitude is very good, which suggests that the one parameter model gluon 
propagator will also provide a good description of hadronic form factors. 

One observes that the mass formula in Eq. ([HP yields an accurate estimate of the mass 
obtained by solving the pion BSE. We find that, with parameters of Eq. fl52|), the right- 
hand-side of Eq. (f4l|) is well described by 

mlNl = 2 (0.45) 3 m R + (2.6) 2 m\ + 150 m 3 R (53) 

on the range m# G [0, 0.02] GeV, from which one may infer a value of (qq) p = —(0.45 GeV) 3 . 
At the value of in Eq. ( |52"D the term linear in ttir contributes almost 96% of the total. 
We see, therefore, that Eq. ([II]) entails m 2 oc m#, for small m^, but that the constant of 
proportionality is not given by the usual definition of the vacuum quark condensate, Eq. ([13]) . 

Our one parameter model for the gluon propagator explicitly excludes the ln[/c 2 ]- 
corrections associated with the anomalous dimensions in QCD. It is therefore inappropriate 
to directly compare m^(/x) in Eq. ( |52|) with the QCD evolution of the commonly quoted 
value of m /J=1GeV ~ 7.5 MeV. j43| (This entails that the same is true of (qq) p -) We note 
that replacing (wd)/k 2 by a^°~ loop \k 2 ) / 'k 2 in Eq. ([I]) would lead to a suppression of the tail 
of the quark mass function, thereby requiring a larger value of ttir to reproduce the pion 
mass and a commensurate change in m t . This represents a quantitative improvement of our 
model but would not change its qualitative features. 

A. p-meson observables. 

We have employed our model gluon propagator in a preliminary study of p-meson prop- 
erties. 

The regularised, generalised ladder approximation to the p-meson BSE is 

F;(p; P) = Z 2 3 C 2 (R) f g 2 D^(p - q) -Ltr [ la i llx S(q+) iT a {P) S{qJ) lv \ F p (q; P) , (54) 

where [T p (P) = 7 M + 7 • PP^/m 2 }. The subtractively renormalised amplitude is given 
by F p (p; P) = Fp(p] P) — F' p (n', P)- We neglected the other Dirac-structures allowed in the 
vector-meson Bethe-Salpeter amplitude. For the p-meson the error introduced by this trun- 
cation is approximately 10%. []36]] The p- and cu-mesons are degenerate at this level of ap- 
proximation. As for the pion, we project this equation onto the lowest Tschebyshev moment 
and solve for Fo(p 2 ,P 2 ), neglecting the coupling to the higher moments. This is a good 
approximation for the p-meson. [2H] 



In this preliminary study we have only solved the quark DSE at real-p 2 . For an on-shell 
p-meson P 2 < and hence Eq. (^) samples the quark propagator at complex values of p 2 . 
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To obtain an approximate solution of Eq. (|54|), without solving the quark DSE at complex- 



p 2 , we introduced an eigenvalue, A(P 2 ), on the right-hand-side of Eq. (54) and solved this 
equation at spacelike values of P 2 , thereby obtaining A(P 2 > 0). For < P 2 < 10 fm~ 2 the 
results could be described by the quadratic (in P 2 ): 

A(P 2 ) = 0.44 - 0.021 P 2 + 0.000076 P 4 (55) 

with a standard-deviation of 0.000044. We compared this with both linear and cubic fitting 
forms: it provides a smaller standard-deviation than the linear form and is monotonic, 
whereas the cubic is not. The value of P 2 for which this algebraic form of A(P 2 ) = 1 
provides the mass estimate presented in Table |. 

The calculated p-meson Bethe-Salpeter amplitude is much narrower in momentum space 
than that of the pion, in agreement with the results of Ref. EH . 



The calculation of g piT7T proceeds in a similar manner. In generalised impulse approx- 
imation the pjiTT coupling can be expressed in terms of a nonlocal coupling functional, 
N p (p,q), which is discussed in Ref. flj?fl . This expression is used to evaluate g P mr{P 2 ) on 
< P 2 < 10 fin -2 . The results were fitted and extrapolated to the calculated mass-shell 
point. The best fit was obtained with: 

g pn7T (P 2 ) = 1.15 - 0.076 P 2 + 0.0013 P 4 - 0.000022P 6 , (56) 

giving a standard- deviation of 0.00023. This form provides a smaller value of the standard- 
deviation than either a linear or quadratic form and is monotonic whereas the quartic is not. 
The value obtained at the calculated on-mass-shell point is given in Table Q. 

These calculations are only a first step. They serve merely to indicate that our one 
parameter model gluon propagator, which was fitted to pion observables, can reasonably be 
expected to provide a good description of other observables too. 



B. Confinement. 

We have also solved the quark DSE for real-p 2 < 0. There is no singularity on the real-p 2 
axis. The solution therefore does not have a Lehmann representation and hence may be 
interpreted as describing a confined particle. 

A plot of l/[p 2 + M(p 2 ) 2 }, which for a free particle would have a pole at the mass-shell 
point, has a broad resonance-like peak centred on p 2 w —(0.55) GeV 2 . This admits an 
interpretation as the "constituent-quark-mass" in our model. 

The form of our solution is suggestive of a pair of complex conjugate poles or branch 
points with timelike real parts and large magnitude imaginary parts. We have made no 
attempt to confirm this. A thorough study must identify whether this structure is an 
artifact of the rainbow approximation, which is known to be associated with unexpected 
behaviour of the fermion propagator in the complex plane |48|-[52| that is modified when the 
vertex is dressed. [|53l 



VI. SUMMARY AND CONCLUSIONS 

Using a confining, one parameter model form for the gluon propagator, Eq. (|4]), which 
incorporates the essence of the solution of realistic, approximate gluon Dyson-Schwinger 
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equations (DSEs), we solved the renormalised, rainbow approximation quark DSE and sub- 
sequently the renormalised, generalised ladder approximation tt- and p-meson Bethe-Salpeter 
equations (BSEs). We varied the parameter in the gluon propagator, m t , which is a mass 
scale that marks the transition between the perturbative and nonperturbative domains, and 
the renormalised current quark mass and obtained a good description of a range of tt- and 
p-meson observables. The value of m t was not known a priori. Good agreement with the 
data required m t ~ 700 MeV, which corresponds to a length of ~ 0.3 fm. 

In studying the pion BSE we were led to a mass formula for the pion, Eq. fl4T|), expressed 
solely in terms of the massive and massless quark propagators. This formula provides a very 
accurate estimate of the pion mass. It is valid to all orders in the renormalised current 
quark mass, and for m# < 20 MeV the nonlinear terms provide a contribution of no more 
than ~10%. 

We obtained numerical solutions of the quark DSE on the timelike-p 2 axis, which showed 
the quark propagator to have no singularity on the real-p 2 axis in our model. We found 
evidence to suggest that, as a function of p 2 , the quark propagator has a pair of complex 
conjugate poles or branch points with timelike real parts and large imaginary parts. Such 
a propagator does not have a Lehmann representation and admits the interpretation of 
describing a confined particle. 

Our study illustrates the manner in which the DSEs can be used to develop a semi- 
phenomenological approach to QCD that incorporates the perturbative, large spacelike-/c 2 
behaviour known from renormalisation group studies and, via an economical parametrisa- 
tion, extrapolates this into the nonperturbative, small spacelike-fc 2 domain. This efficacious, 
nonperturbative approach allows for the correlation of a large range of observables via very 
few parameters, which it may be possible to relate to the fundamental parameters of QCD. 
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TABLES 





Calculated 


Experiment 


^Mass— Formula: J2Jq] 


138.7 MeV 


138.3± 0.5 


™.Mass— Formula: f 77 1 


137.2 




^ BS— Equation 


139.5 




fn\F ] 


92.4 MeV 


92.4± 0.3 


fn\B } 


92.3 






102 




/I |_ U J /I \_ UJ 


0.24 


0.31± 0.004 


IK U J 


0.16 


0.21± 0.02 




-0.041 


-0.040 ± 0.003 


«M 


0.028 


0.038 ± 0.003 


a°[F ] 


0.0022 


0.0017 ± 0.0003 




0.0013 




flV°77 [-^o] 


0.45 


0.50 ± 0.02 




0.971 GeV 


0.770 ± 0.001 


^7 p7T7T [Fq] 


4.07 


6.07 ± 0.02 



TABLE I. Observables calculated using the parameter values in Eq. (p2j). The experimental 
values of the tt-tt scattering lengths are discussed in Refs. pjfL^I- The other experimental values are 
are taken from Ref. [ )43| | . [Bq] indicates that the quantity was calculated using the approximation 
of Eq. (^7j) while [Fq] indicates it was calculated using the zeroth order Tschebyshev moment 

is discussed in 



obtained in a direct solution of the BSE, Sec. IV A, The anomalous coupling g n o^ 
Ref. ||. See Sec. V A for a discussion of the p-meson observables. The difference between N n and 
/tt is a measure of the accuracy of the approximation of Eq. (^7|). That between the calculated 



and experimental values of r n N^ is a measure of the importance of final-state tt-tt interactions and 
photon- /3-meson mixing. [44] Final-state tt-tt interactions are also neglected in the ca 
the scattering lengths [Q] and g pn7T . Pion-loop corrections to m p are of the order of 5% 
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